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ABSTRACT 

We present the discovery of two transiting exoplanets. HAT-P-28b orbits a V=13.03 G3 dwarf star 
with a period P = 3.2572 d and has a mass of 0.63 ± 0.04 Mj and a radius of i?j yielding 

a mean density of 0.44 ± 0.09 gcm~ 3 . HAT-P-29b orbits a V=11.90 F8 dwarf star with a period 
P = 5.7232 d and has a mass of 0.78^04 -^J an d a radius of l.ll^Jg-Rj yielding a mean density of 
0.71 ±0.18 gem -3 . We discuss the properties of these planets in the context of other known transiting 
planets. 

Subject headings: planetary systems — stars: individual (HAT-P-28, GSC 2284-00503, HAT-P-29, 
GSC 3293-01539) techniques: spectroscopic, photometric 



1. INTRODUCTION 

When a planet transits its host star, we can detect 
not only the presence of the planet, but also learn a 
great deal about its physical properties: we can infer the 
planet's mass and radius and thus its bulk density and 
com position, its temperature a nd atmospheric composi- 
tion (jCharbonneau et al.ll2005l ) , and the dynamics of the 
planet ary system such a s the projected spin-orbit align- 
ment (W mn eTaDHool). With over 100 transiting exo- 
planets (TEPs) discovered in wide-angle transit surveys 
and the wealth of transiting planet candidates from the 
Kepler mission ([Borucki et al.lf20TlT ) . statistical analyzes 
are beginning to support constraints on the planet for- 
matio n models (jTriaud et al.l I2010L iMorton fc Johnsonl 
l20ll . 

Here we report the discovery of two new TEPs, HAT-P- 
28b and HAT-P-29b, orbiting the stars GSC 2284-00503 
and GSC 3293-01539, respectively. The Hun garian-made 
Automated Telescope Network (HATNet; IBakos et al.1 
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2004) has been in operation since 2003, and has now cov- 
ered approximately 14% of the sky, searching for TEPs 
around bright stars (9.5 < r < 14.5). Although HAT- 
P-28b and HAT-P-29b are of the most common types 
of confirmed TEPs with respect to their mass and radii, 
they contribute to the characterization of the growing 
population of Jupiter-sized planets known to transit. 
The two discovered planets reported here orbit relatively 
faint stars compared to the current population of vali- 
dated transiting planets, but the typical HATNet planets 
have relatively bright host stars compared to the wealth 
of planets candidates discovered by Kepler, and are hence 
more amenable to follow-up confirmation and character- 
ization. 

In Section [2] we briefly describe the initial photometric 
detection and subsequent spectroscopic and photomet- 
ric follow-up observations used to confirm the planetary 
nature of HAT-P-28b and HAT-P-29b. In Section [5] we 
describe the data analysis which we performed to deter- 
mine the planetary and stellar parameters. Our findings 
are discussed in Section |4] 

2. OBSERVATIONS 

HAT-P-28b and HAT-P-29b were discovered by HAT- 
Net following an observational procedure which has been 
descr ibed in detail in several previous discove ry papers 
(e.g. IBakos et al.l l2010t iBuchhave etafl [2010l ) . Briefly, 
the procedure consists of the following steps: 

1. Photometric observations with one or more of the 
HATNet telescopes to detect stars with periodic 
transit-like light curves. 

2. High-resolution, low-S/N "reconnaissance" spec- 
troscopy observations to reject many false positives 
(such as M dwarf stars transiting G or F stars, or 
eclipsing binary stars blended with brighter giant 
stars). 

3. Higher-precision photometric observations during 
transit to confirm the detection, show that it is 
consistent with a transiting planet, and use in de- 
termining the system parameters. 



2 



Buchhave et al. 



4. High-resolution, high-S/N "confirmation" spec- 
troscopy to detect the orbital motion of the star 
due to the planet, characterize the host star, and 
rule-out subtle blend scenarios. 

Below we describe details of this procedure that are per- 
tinent specifically to the discovery of HAT-P-28b and 
HAT-P-29b. 

2.1. Photometric detection 
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Fig. 1.— Unbinned light curve of HAT-P-28 including all 10,400 
instrumental i?-band 5.5 minute cadence measurements obtained 
with the HAT-6, HAT-7, HAT- 8 and HAT-9 telescopes of HAT- 
Net (see Table ft}, and folded with the period P = 3.257215 ± 
0.000007 days resulting from the global fit described in Section \3\ 
The soli d lin e shows a simplified transit model fit to the light curve 
(Section l3.2H . The lower panel shows a zoomed-in view of the tran- 
sit; the dark filled points show the light curve binned in phase using 
a bin-size of 0.002. 



These observations revealed box-like transit signals in 
the light curves of the stars with an apparent depth of 
~ 16.1 mmag and ~ 9.1mmag, and a period of P = 
3.2572 days and P = 5.7232 days, respectively (see Fig- 
ures [T] and [2]) . 

2.2. Reconnaissance Spectroscopy 

High-resolution, low-S/N reconnaissance spectra were 
obtained for HAT-P-28 and HAT-P-29 using the 
Harvard-Smithsonian C enter for Astro physics (CfA) Dig- 
ital Speedometer (DS; iLathaml I1992T ) and the Tilling- 
hast Reflector Echcllc Spectrograph (TRES; H urcsz 
120081 ) on the FLWO 1.5 m telescope as well as the 
Fiber-fed Echelle Spectrograph (FIES) on the 2.5 m 
Nordic Optical Telescope ( NOT) at La Palma, Spain 
(jDiupvik fc Andersen! [2010h . The reconnaissance ob- 
servations are listed in Table [2] and were used to esti- 
mate the effective temperature, surface gravity, and pro- 
jected rotational ve locity of the host stars as described by 
iTorres et al.l ([2003)- The procedures used for reduction 
and extraction of high precisi on radial velocit i es of the 
FIES spectra are described in lBuchhave et afl (|2010h . 

The CfA DS reconnaissance observations revealed no 
detectable RV variation at the ~ lkms -1 precision of 
the observations. The TRES and FIES reconnaissance 
observations showed radial velocities with a low RMS 
(HAT-P-28b showed an RMS of 12 ms" 1 and the HAT- 
P-29 showed an RMS of 82 ms _1 ), which both could be 
consistent with a planetary companion. Additionally the 
spectra are consistent with the host stars being single, 
slowly-rotating, dwarf stars. 
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Fig. 2.— Unbinned light curve of HAT-P-29 including all 3,453 
instrumental Sloan r 5.5 minute cadence measurements obtained 
with the HAT-5 and HAT-8 telescopes and folded with the period 
P = 5.723186 ± 0.000049 days. See Figure [TJ caption for details. 

Table [T] summarizes the photometric observations 
of HAT-P-28b and HAT-P-29b. Altogether 10,400 
instrumental i?-band 5.5 minute cadence photomet- 
ric measurements of HAT-P-28 (GSC 2284-00503; 
2MASS 00520018+3443422; a = 00 h 52 m 00.18s 6 = 
+34 43'42.3"; J2000; V=13.03. lDroege et al. 2006) were 
obtained with the HAT-6, HAT-7, HAT-8 and HAT- 
9 telescopes. 3,453 Sloan r measurements of HAT-P- 
29 (GSC 3293-01539; 2MASS 02123147+5146435; a = 
02 h 12 m 31.46s, 5 = +51°46'43.5"; J2000; V=11.90) were 
obtained with the HAT-5 and HAT-8 telescopes. 



2.3. Photometric follow-up observations 

In-transit photometric observations of HAT-P-28 and 
HAT-P-29 were obtained with the KeplerCam CCD cam- 
era on the FLWO 1.2 m telescope and with the 2.0 m 
Faulkes Telescope North (FTN) at Haleakala Observa- 
tory in Hawaii. These observations are summarized in 
Table [J 

The images were reduced to light curves following the 
apertu re photometry procedure described bv lBakos et alJ 
(2010). Fi gures [3] and [4] show the final trend-filtered 
light curves together with the best-fit transit model (see 
Section l3"T2"j) ; the individual measurements are provided 
in Table [3] and H 

2.4. High resolution, high S/N spectroscopy 

We obtained high-resolution, high-S/N spectra of 
HAT-P-28 using NOT/FIES as well as HIRES 
(|Vogt et alJ 119941) mounted on the Keck I telescope on 
Mauna Kea, Hawaii. A total of 10 high S/N FIES obser- 
vations were obtained between 3 and 18 October 2009, 
and a total of 8 HIRES observations, including 6 through 
an iodine gas absorption cell and two without the cell, 
were obtained between 1 October 2009 and 23 February 
2010. For HAT-P-29, 9 HIRES spectra were obtained be- 
tween 26 September and 14 December 2010, one of which 
was without the iodine cell. 

For the FIES observations we used the medium- 
resolution fiber, giving a resolution of A/AA « 46,000 
and a wavelength coverage of -3600-7400 A. The spec- 
tra were bracketed by ThAr calibration frames used to 
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TABLE 1 

Summary of photometric observations 



Instrument /Field a 



Date(s) 



Number of Images Cadence (s) Filter 



HAT-P-28b 

HAT-7/163 

HAT-8/163 

HAT-6/164 

HAT-9/164 

KeplerCam b 

KeplerCam b 

FTN C 

HAT-P-29b 

HAT-5/089 
HAT-8/089 
KeplerCam b 
KeplerCam b 



2007 Sep-2008 Jan 
2007 Sep-2008 Jan 
2007 Sep-2008 Feb 
2007 Sep-2008 Feb 
2010 Sep 03 
2010 Sep 06 
2010 Oct 02 



2008 Oct-2009 Mar 
2008 Oct-2009 Mar 
2011 Jan 01 
2011 Jan 24 



2329 
1624 
3685 
2779 
174 
129 
476 



2779 
674 
199 
261 



330 
330 
330 
330 
133 
163 
35 



330 
330 
73 
43 



R 
R 
R 
R 

Sloan i 
Sloan i 
Sloan i 



Sloan r 
Sloan r 
Sloan i 
Sloan i 



HAT-P-29 is in one 



a HAT-P-28 is in two HATNet fields, internally labeled as 163 and 164 
field, labeled 089. 

b Observations from the KeplerCam CCD camera on the FLWO 1.2 m telescope 

c Observations from the 2.0 m Faulkes Telescope North (FTN) at Halcakala Observatory in 

Hawaii 



TABLE 2 

Summary of reconnaissance spectroscopy observations 



Instrument Date 


T e ff* 


logs* 


v sin i 


7RV a 


(K) 


(cgs) 


(kms -1 ) 


(km s~ - 1 ) 




HAT-P-28b 










DS 13 Jan 2009 


5750 


1.5 





+44.32 


DS 14 Jan 2009 


5750 


4.5 





+42.34 


FIES 19 Jan 2009 


5750 


4.5 


2 


+43.85 


FIES 21 Jan 2009 


5750 


4.5 


2 


+43.87 


HAT-P-29b 










FIES 17 Aug 2010 


6000 


4.0 


6 


-21.67 


FIES 17 Aug 2010 


5750 


4.0 


6 


-21.76 


TRES 21 Jan 2009 


5750 


4.5 


6 


-21.59 


a The heliocentric RV of the target. 



TABLE 3 

High-precision differential photometry of HAT-P-28 



TABLE 4 

High-precision differential photometry of HAT-P-29 



BJD 

(2,400,000+) 


Mag a 


""Mag 


Mag(orig) b Filter 


BJD 

(2,400,000+) 


Mag a 


"Mag 


Mag(orig) b Filter 


55443.68956 


0.01257 


0.00104 


11.6842 i 


55563.67299 


0.00148 


0.00114 


10.8583 i 


55443.69111 


0.01476 


0.00102 


11.6858 i 


55563.67366 


0.00031 


0.00115 


10.8581 i 


55443.69266 


0.01656 


0.00102 


11.6878 i 


55563.67444 


-0.00043 


0.00099 


10.8565 i 


55443.69421 


0.01308 


0.00102 


11.6852 i 


55563.67528 


-0.00158 


0.00099 


10.8556 i 


55443.69576 


0.01364 


0.00102 


11.6776 i 


55563.67614 


-0.00072 


0.00097 


10.8564 i 


55443.70735 


0.01317 


0.00100 


11.6841 i 


55563.67700 


0.00009 


0.00096 


10.8567 i 


55443.70889 


0.01369 


0.00100 


11.6852 i 


55563.67784 


-0.00134 


0.00094 


10.8549 i 


55443.71044 


0.00874 


0.00099 


11.6816 i 


55563.67870 


0.00242 


0.00095 


10.8600 i 


55443.71197 


0.00762 


0.00100 


11.6786 i 


55563.67955 


0.00023 


0.00094 


10.8577 i 


55443.71663 


0.00139 


0.00100 


11.6848 i 


55563.68042 


0.00100 


0.00095 


10.8581 i 



Note. — This table is available in a machine-readable form in 
the online journal. A portion is shown here for guidance regarding 
its form and content. 

a The out-of-transit level has been subtracted. These magni- 
tudes have been subjected to the External Parameter Dccorrclation 
(EPD) technique and Trend Filtering Algorithm (TFA) procedures 
(sec Bakos ct al. 2010), carried out simultaneously with the transit 
fit. 

Raw magnitude values without application of the EPD and TFA 
procedures. 



Note. — The notes of this table arc identical to the notes in Table [3l 

determine the fiducial wavelength calibration and were 
reduced to relative RVs in the solar system barycentric 
fra me following a cross-c orrelation procedure described 
bv lBuchhave et all (|2Q1(1) . For HIRES, we set the spec- 
trometer slit to 0'/86, resulting in a resolving power of 
A/AA ~ 55,000 with a wavelength coverage of ~3800- 
8000 A. We reduced these spectra to relative RVs in 
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Fig. 3. — Unbinned instrumental Sloan i-band transit light 
curves of HAT-P-28, acquired with KeplerCam at the FLWO 1.2 m 
telescope on 2010 Sep 03 and 2010 Sep 06, and with the Faulkes 
Telescope North on 2010 Oct 02, from top to bottom. The light 
curves have be en c lean ed of instrumenta l variations using the 
ELTG model (§ EHl and [Bakos et afll20T(ih . Curves are displaced 
vertically for clari ty. Our best fit from the global modeling de- 
scribed in Section 13.21 is shown by the solid lines. Residuals from 
the fits are displayed at the bottom, in the same order as the top 
curves. The error bars represent the photon and background shot 
noise, plus the readout noise. 

the solar system barycentri c frame following; a procedure 
based on that described bv lButler etall (|1996f l. The RV 
measurements and their uncertainties are listed in Ta- 
ble [5] and [51 The period- folded data, along with our best 
fit described below in Section [31 are displayed in Figures 

E] and ED 

Following Pforres et all (|2007f ) and lQueloz et al.l(|2001l ). 
w e conducted a bisect or span analysis as described in §5 
of IBakos et all (|2007f) of the Keck and FIES spectra to 
rule out the possibility that the objects are eclipsing bi- 
nary star systems blended with light from a third star 
rather than a planet. As seen in Figures [5] and [6] the bi- 
sector spans show no variation in phase with the orbital 
period and have a total scatter significantly less than 
the RV semi-amplitude, which supports the interpreta- 
tion that the velocity variations are due to a planetary 
companions. 

In the same figures, we also show the calibrated S 
indices, which are a measure of the chromospheric ac- 
tivity of the stars derived from the flux in the cores of 
the Ca II H and K lines (iVauehan. Preston fe Wilsonl 
[197& llsaacson fe Fischerl l2010)." We also measured the 
R' HK indices for the systems to be logi?^^ = —4.984 
and fogi?' w ^ = — 5.105, respectively, as described by 
iNoves et all (|1984[ ). We do not detect any significant 
variation of the indices correlated with orbital phases. 







-0.15 -0.1 -0.05 0.05 0.1 0.15 

Time from transit center (days) 

Fig. 4. — Unbinned instrumental Sloan i-band transit light 
curves of HAT-P-29, acquired with KeplerCam at the FLWO 1.2 m 
telescope on 2011 Jan 01 and Jan 24. See caption in Figurc|3]for 
details. 

3. ANALYSIS 

The analysis of the HAT-P-28 and HAT-P-29 systems, 
including determinations of the properties of the host 
stars and planets, was carried out i n a fashion simila r 
to previous HATNet discoveries (e.g. Bakos e t al.l[2010[) . 
Below we briefly summarize the procedure and results. 

3.1. Properties of the parent star 

Stellar atmospheric parameters were measured using 
our template spectrum obtained with the Keck/HIRES 
instrument, and the analysi s package known as Spec 
troscopy Made Easy (SME; IValenti fe Piskunovl 119961 
along with the atomic line database of IValenti fc Fische 
(2005). SME yielded initial values and uncertainties 
(which, for the first two parameters, we have conserva- 
tively doubled to include our estimates of th e systematic 
error s based on prior experience (see e.g. INoves et al.l 
120081 ) ). For HAT-P-28 we found the following param- 
eters: effective temperature T e a* — 5680 ± 90 K, metal- 
licity [Fe/H] = +0.12 ± 0.08 dex, projected rotational 
velocity wsini = 0.2+Q2 kms^ 1 , and stellar surface grav- 
ity log 3* = 4.42 ± 0.06 (cgs). For HAT-P-29 we found: 
effective temperature T c g+ — 6087 ± 88 K, metallicity 
[Fe/H] = 0.21 ± 0.08 dex, projected rotational velocity 
vsini = 3.9 ± 0.5kms _1 , and stellar surface gravity 
log = 4.31 ±0.06 (cgs). 

Our initial values of T c ff*, logg*, and [Fe/H] were 
used to determine the limb-darkening coefficients needed 
in the global modeling of the follow-up photometry. 
This modeling is briefly summarized in Sec tion 13.21 with 
more details given by B akos et al.l ((20101 ) . Following 
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TABLE 5 

Relative radial velocities, bisector spans, and activity index measurements 

of HAT-P-28. 



BJD 


RV a 




BS 


OF3S 




s c 


Phase 


Instrument 


(2,400,000+) 


(m s" 1 ) 


(ms ) 


(ms- 1 ) 


(ms ) 












55107.11869 


83.86 


3.76 


—8.78 


3.78 





177 





680 


Keck 


rri no 171 1 m 

55108.71101 


—72.31 


12.10 


18.20 


13.60 









168 


FIES 


5511U. 08494 


77.09 


12.60 


— 11.60 


17.90 









774 


FIES 


ooiii.oyy ( i 


—22.21 


9.10 


—3.80 


11.40 









055 


1 lljb 


55112. 63861 


—57.41 


12.60 


9.10 


12.30 









374 


FIES 


55110.00539 


78.79 


8.70 


—9.40 


7.70 









671 


FIES 


55114.50809 


10.09 


13.90 


10.90 


9.70 









967 


1 lljb 


55115.56626 


-99.31 


10.80 


6.00 


10.70 









273 


FIES 


55116.55720 


34.79 


10.00 


1.10 


13.80 









577 


FIES 


55122.55573 


-43.01 


8.70 


-14.90 


12.60 









419 


FIES 


55123.50498 


61.29 


11.90 


1.40 


11.30 









710 


FIES 


55125.49641 


-68.71 


9.70 


-6.90 


7.50 









322 


FIES 


55191.86077 






2.33 


2.14 





163 





696 


Keck 


55191.87545 


80.26 


2.74 


8.43 


2.12 





168 





701 


Keck 


55193.81948 


-83.44 


2.85 


10.12 


3.37 





168 





298 


Keck 


55198.83008 


64.79 


2.92 


7.03 


2.85 





166 





836 


Keck 


55229.73714 


-99.43 


3.16 


-3.11 


3.96 





158 





325 


Keck 


55251.74043 






-0.42 


3.40 





154 





080 


Keck 


55251.75531 


-64.03 


4.85 


-15.61 


11.55 





102 





085 


Keck 



Note. — For the iodinc-frcc template exposures there is no RV measurement, but the BS 
and S index can still be determined. 

a The zero-point of these velocities is arbitrary. An overall offset 7 rG i fitted separately to the 
FIES and Keck velocities in Scction |3.2| has been subtracted. 

Internal errors excluding the component of astrophysical/instrumental jitter considered in 
Section [3T2] 

c Relative chromospheric activity index, calibrated to the scale of 

IVaughan. Preston fe Wilson! J 1978ft . 



TABLE 6 

Relative radial velocities, bisector spans, and activity index measurements 

of HAT-P-29. 



BJD 


RV a 


_ b 
CTRV 


BS 


"FJS 


S c 


Phase 


Instrument 


(2,400,000+) 


(ms — 1 ) 


(ms — 1 ) 


(ms" 1 ) 


(ms" 1 ) 








55465.99028 






-0.20 


1.44 


0.148 


0.900 


Keck 


55466.00284 


37.80 


2.16 


-3.72 


2.05 


0.145 


0.902 


Keck 


55467.94293 


-85.06 


2.47 


-3.19 


1.42 


0.142 


0.241 


Keck 


55469.02138 


-15.21 


2.31 


2.73 


1.87 


0.140 


0.429 


Keck 


55469.92851 


42.10 


2.53 


-3.47 


3.07 


0.140 


0.588 


Keck 


55490.91798 


-79.20 


2.51 


-9.19 


1.97 


0.143 


0.255 


Keck 


55500.96329 






3.48 


1.14 


0.133 


0.010 


Keck 


55521.91781 


80.29 


2.39 


-6.04 


7.56 


0.147 


0.672 


Keck 


55523.00463 


45.97 


2.49 


10.14 


5.01 


0.147 


0.862 


Keck 


55544.95590 


68.78 


3.04 


-7.18 


3.45 


0.136 


0.697 


Keck 



NOTE. — For the iodinc-frcc template exposures there is no RV measurement, but the BS 
and S index can still be determined. 

a The zero-point of these velocities is arbitrary. An overall offset 7 rG i fitted separately to the 
FIES and Keck velocities in Section |3.2| has been subtracted. 

k Internal errors excluding the component of astrophysical/instrumental jitter considered in 
Section [3]2] 

c Relative chromospheric activity 

IVaughan. Preston & Wilsonl d 1 978ft . 

iSozzetti et aLl(|2007l ) we used a / R+ (the normalized semi- 
major axis), Tcff* and [Fe/H] togethe r with the Yon sei- 
Yale (YY) stellar evolution models (|Yi et al.l l200l to 
determine probability distributions of other stellar prop- 
erties (such as logg+, and M*). Note that a/R+ 
is related to the mean stellar density, and is determined 
from the analysis of the light curves and RV curves. This 
procedu re has been described in further detail by iPall 
(|2009bD . 

The inferred location of the stars i n a di agram of a/i?* 
versus T e fj* together with lYi et al.l (|2001[ ) isochrones, is 



index, calibrated to the scale of 

shown in Figures [7] and The SME results for the sur- 
face gravity yielded log 5* = 4.36 ± 0.06 for HAT-P-28 
and log = 4.34 ± 0.06 for HAT-P-29. Both values are 
consistent with our initial SME analysis and we therefore 
did not carry out a second iteration of SME and adopted 
these values for log g+ and the values of T e ff+, [Fe/H], and 
v sin i stated above as the final atmospheric properties of 
the star. The stellar parameters are listed in Table[7l We 
find that HAT-P-28 has a mass M* = 1.025 ± 0.047 M , 
radius R„ = 1.103±g-g|| R e , and an age of 6.1±f-| Gyr 
and HAT-P-29 has a mass M* = 1.207 ± 0.046 Af Q , ra- 
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Fig. 5. — Top panel: High-precision relative RV measurements 
for HAT-P-28 from Keck/HIRES (filled circles) and NOT/FIES 
(open triangles) shown as a function of orbital phase, along with 
our best-fit model (see Table [SJl . Repeated points are shown as 
open circles and crosses for Keck/HIRES and NOT/FIES, respec- 
tively. Zero phase corresponds to the time of mid- transit. The 
center-of-mass velocity has been subtracted. Second panel: Ve- 
locity O — C residuals from the best fit. The error bars of the 
Keck/HIRES measurements include a component from astrophysi- 
cal/instrumcntal jitter (15.6 ms ; see Section l3,2l l; it was unnec- 
essary to add jitter to the NOT/FIES measurements. The RMS 
of the residuals is 16.1ms" 1 and 10.1ms -1 for the Keck/HIRES 
and NOT/FIES observations respectively. Third panel: Bisector 
spans (BS), with the mean value subtracted. The measurement 
from the template spectrum is included. Bottom panel: Relative 
chromospheric activity index S measured from the Keck spectra. 
Note the different vertical scales of the panels. 

dius R* = 1.224±° ^i? , and an age of 2.2 ± 1.0 Gyr, 
according to these models. 

As a check on the stellar evolution modeling we com- 
pare the observed photometric color indices of the stars 
to the values predicted by the modeling. We use 
the ne ar-infrared magnitud es from the 2MASS Cata- 
logue (jSkrutskie et all 120061 ). which for HAT-P-28 are 
J 2 mass = H-561 ± 0.020, H 2MAS s = 11.206 ± 0.023 
and -£T 2 mass = 11.104 ± 0.018 and for HAT-P-29 are 
J 2 mass = 10.648 ± 0.023, i? 2 MASS = 10.395 ± 0.023 and 
#2MASS = 10.297±0.020. Converting these to the photo- 
metric s ystem of th e models (ESO) using the transforma- 
tions bv lCarpenterl ()2001l ) we find a measured color index 
of J - K = 0.487 ± 0.030 and J — K = 0.375 ± 0.033, 
respectively. This is within 3cr of the value predicted 
from the isochrones of J — K = 0.410 ± 0.020 and 
J-K = 0.340 ± 0.020, respectively. 

From the observed and predicted J — K we find 
E(J-K) = 0.077± 0.036 and E{J-K) = 0.035 ±0.039, 
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Fig. 6. — Top panel: High-precision relative RV measurements 
for HAT-P-29 from Keck/HIRES shown as a function of orbital 
phase, along with our best-fit model (see Table [8|. Zero phase 
corresponds to the time of mid-transit. The center-of-mass velocity 
has been subtracted. Repeated points are shown as open circles. 
Second panel: Velocity O—C residuals from the best fit. The error 
bars of the Keck/HIRES measurements include a component from 
astrophysical/instrumental jitter (6.0ms -1 ; see Section [32} . The 
RMS of the residuals is 6.5ms" 1 . Third panel: Bisector spans 
(BS), with the mean value subtracted. Bottom panel: Relative 
chromospheric activity index S measured from the Keck spectra. 
Note the different vertical scales of the panels. 
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Fig. 7. — Model isochrones from Yi et al. (2001) for the measured 
metallicity of HAT-P-28, [Fe/H]= +0.12, and ages between 1 and 
14 Gyr in 1 Gyr increments (left to right). The adopted values of 
T e ff* and a/-R* are shown together with their 1<t and 2<t confidence 
ellipsoids. 



respectively, and using ICardelli et all (fl989h we find 
A{K) = 0.052 ± 0.024 and A{K) = 0.024 ± 0.026, re- 
spectively, which gives a reddening corrected distance 
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Effective temperature [K] 

Fig. 8. — Model isochrones from Yi ct al. (2001) for the measured 
metallicity of HAT-P-29, [Fe/H]= 0.21, and ages between 1 and 
13Gyr in 1 Gyr increments (left to right). The adopted values of 
Tefj* and a/ R+ are shown together with their 1<t and 2<r confidence 
ellipsoids. 

estimate, using the K magnitude, of 395^26 P c f° r HAT- 
P-28 and 322^ pc for HAT-P-29. The uncertainty does 
not include difficult to quantify systematics in the model 
isochrones. 

3.2. Global modeling of the data 

We modeled the HATNet photometry, the follow-up 
photometry, and the high-precision RV m easurements 
using t he p rocedure described in detail by iBakos et al.l 
(201(3); El ([2009). To summarize, we s i multa neously 
fit these data using the iMandel fc Agoll (|2002t) model 
for the follow-up light curves from KeplerCam and 
FTN, quadratic li mb dark e ning c oefficients interpolated 
from the tables by G]arct (2004), a simplified no- limb- 
darkening transit model for the HATNet light curves, 
and a Keplerian orbit for the RVs. 

The free parameters in the fit include the time of the 
first and last transit center observed, the normalized 
planetary radius p = the square of the impact 

parameter b 2 , the reciprocal of the half duration of the 
transit C/ii*, the RV semiamplitudc K, Lagrangian el- 
ements k = ecosu, and h = esinw, the HATNet blend 
factors (which account for possible dilution of the transit 
in the HATNet light curves from background stars due 
to the broad PSF (~ 26" FWHM) and the application 
of trend-filtering without signal reconstruction) , the out- 
of-transit magnitudes for each HATNet field, the relative 
zero-points 7 rc i, hires and 7rei,FiES of the HIRES/Keck 
and FIES/NOT RVs, and 22 instrumental parameters 
from the "ELTG" model for systematics. 

We inflated the errors for the HIRES velocities by 
adding in quadrature a jitter of 15.6 ms -1 for HAT- 
P-28b and 6.0 ms" 1 for HAT-P-29b to the formal er- 
rors to achieve a x 2 /dof = 1. The internal uncertainty 
estimate of the FIES velocities, which is estimated by 
a = RMS(v)/y/N, where v is the RVs of the individual 
order s and N is the number of orders (see lBuchhave et al.l 
l2010f ) , might tend to slightly overestimate the uncertain- 
ties of the individual measurements and thus we do not 
need to add jitter to the FIES velocities to achieve a 
X 2 /dof=l. 

The derived planetary parameters can be found in Ta- 
ble 1 We find: 

• HAT-P-28b - The planet has a mass of M p = 



0.626 ± 0.037 Mj and a radius of R p = 
1.212lo;Jg| Rj, a mean density of p p = 0.44 ± 
0.09 gem" 3 . 



• HAT-P-29b 



The planet has a mass of 



M p = 0.778±g;^o Mj and a radius of R p = 



l.lO?!";^ a mean density of p p 
0.18gcm- 3 . 



0.71 ± 



We find an eccentricity of e = 0.051 ± 0.033 and 
e = 0.095 ± 0.047 for HAT-P-28b and HAT-P-29b, re- 
spectively, which are both consistent with a circular or- 
bit. 

4. DISCUSSION 

We present the discovery of the transiting planets 
HAT-P-28b and HAT-P-29b. Figure M shows the two 
planets on a mass-radius diagram comparing them to 
the other known TEPs. We discuss the properties of 
each planet below. 

4.1. HAT-P-28b 

HAT-P-28b has a period of P = 3.257215 ±0.000007d, 
a mass of M p = 0.63 ± 0.04 Mj and a radius of R p = 
1.21±Q;Jg leading to a mean density p p — 0.44 ± 
0.09 gem -3 . We find an eccentricity of e = 0.051 ±0.033, 
which is consistent w ith a circular orbit. Both Kepler-6b 
(jDunham et al|2010D with M v = 0-67 Mj, R p = 1.32 
P = 3.23 d Kipping fc Bakosl [20Tol) and Kepler-8b 



(| Jenkins et al.ll2010D with M p = 0.6 Mj, R p = 1.42 



P = 3.52 d (|Kipping k, Bakosll2010[ ) are quite similar to 
HAT-P-28b in terms of their mass and periods, but have 
slightly larger radii. A significant number of other tran- 
siting planets have very comparable properties to HAT- 
P-28b. In fact, in histograms of mass and radius of the 
discovered transiting planets with a bin-size of 0.5 Mj 
and 0.25 Rj, respectively, HAT-P-28b resides in the most 
populated bins. 

We h ave compared HAT-P -28b to the theoretical mod- 
els from lFortnev et ail (|2007| ) by interpolating the models 
to the solar equivalent semimajor axis of a = 0.0434 ± 
0.0007 AU, the result of which can be seen overplotted 
in Figure [9l Given the estimated age of 6.ll 2 /g Gyr, we 
find that HAT-P-28b may be slightly larger than what 
is allowed by the theoretical models, but otherwise con- 
sistent with models of a H/He-dominated planet with a 
negligible core mass. 

4.2. HAT-P-29b 
HAT-P-29b has a period of P = 5.723186±0.000049d, 



a mass of M p 

1+0.14 



0.78 



-0.08 
-0.04 



Mj and a radius of R p 



1.11_ 08 leading to a mean density p p — 0.71 ± 

0.18 gem -3 , and has an eccentricity of e = 0.095 ±0.047, 
which is consistent with a circular orbit. Again, HAT- 
P-29b is quite similar to a number of other transiting 
planets and is almost identical in mass and radius to 
TrES-1 (Mr, = 75 Mj, R p = 1.08 P = 3.03 d, 
lAlonso et al.l 120041 ) . albeit with a slightly longer period. 

When compared with the iFortney et al.l (|2007f ) plane- 
tary models, HAT-P-29b could be consistent with a 4.5 
Gyr model for a H/He-dominated planet with a negligi- 
ble core mass but also with a 1 Gyr model with a core 
mass of M c = 10 M m . 
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TABLE 7 

Stellar parameters for HAT-P-28 and HAT-P-29 



Parameter HAT-P-28 

Spectroscopic properties 

T eff * (K) 

[Fe/H] 

vsini (kms 1 ) . . . 

"mac (kms -1 ) 

i) mic (kms -1 ) 

7rv (kms -1 ) 

Photometric properties 

V (mag) 

V-I c (mag) 

J (mag) 

H (mag) 

K s (mag) 

Derived properties 

M* (M ) 

H* (Re) 

log g* (cgs) 

L* (Lq) 

M v (mag) 

M K (mag.ESO) . . 

Age (Gyr) 

Distance (pc) 395 



HAT-P-29 



Source 



5680 ± 90 


6087 ± 88 


SME a 


+0.12 ±0.08 


0.21 + 0.08 


SME 


2+ - 5 

u ' z -0.2 


3.9 + 0.5 


SME 


3.88 


4.50 


SME 


0.85 


0.85 


SME 


+43.884 ± 0.068 


-21.670 + 0.08 


FIES 


13.03 


11.90 


TASS 


0.85 + 0.18 


0.85 ± 0.18 


TASS 



11.561 + 0.020 
11.206 + 0.023 
11.104 + 0.018 



1.025 + 0.047 
1 103+ 091 
4.36 + 0.06 

i ,,+0.23 

4.70 + 0.19 
3.11+0.16 

T+2.6 
-1.9 



10.648 ± 0.023 
10.395 ± 0.023 
10.297 ± 0.020 



1.224" 1 



6.r 



1.207 + 0.046 

4+0.133 
-0.075 

4.34 + 0.06 

4.11 ± 0.21 
2.77 + 0.19 
2.2 ± 1.0 



+34 
-26 



322 



+35 



2MASS 
2MASS 
2MASS 



YY+a/i?* 
YY+a/i?* 
YY+a/fl* 
YY+a/iJ* 
YY +a/R„ 
YY '+a/R t 
YY +a/R t 
YY +a/R+ 



hSME b 

r-SME 

hSME 

rSME 

S-SME 

rSME 
r-SME 

r-SME 



a SME — "Spectroscopy Made Easy" package for the analysis of high-resolution spectra 
iValentP fe Piskunov 199(|). These parameters rely primarily on SME, but have a small 
dependence also on the iterative analysis incorporating the isochrone search and global 
modeling of the data, as described in the text. 

b YY+o/fl,-fSME = Based on the YY isochrones iYi et al.ll2001f) . a/R* as a luminosity 
indicator, and the SME results. 
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TABLE 8 
Orbital and planetary parameters 



Parameter HAT-P-28b HAT-P-29b 

Light curve parameters 

P (days) 3.257215 ±0.000007 5.723186 ± 0.000049 

T c (BJD) a 2455417.59832 ± 0.00053 2455197.57540 ± 0.00181 

T14 (days) a 0.1341 ± 0.0020 0.1407 ± 0.0074 

T12 = T 34 (days) a 0.0148 ± 0.0015 0.0177 ± 0.0024 

a/R t 8.43 ±0.57 11.70^ ^ 

C/R-k 16.77 ±0.16 16.22 ±0.93 

R p /R* 0.1130 ±0.0024 0.0927 ± 0.0028 

0.089T°;°g 0-34S±g : g5g 

6 = acosi/-R* 0.299l°;j23 0.59ll°;°g4 

i (deg) 88.0 ±0.9 87.1+qJ 

Limb-darkening coefficients b 

Oj (linear term) 0.2879 0.2273 

bi (quadratic term) 0.3220 0.3581 

RV parameters 

K (ms" 1 ) 84.7 ±4.2 78.3 ± 5.9 

ecosaj c -0.010 ±0.028 -0.084^° °^ 

esinu; c -0.022 ± 0.053 0.016 ± 0.058 

e 0.051 ±0.033 0.095 ± 0.047 

w (deg) 233 ± 90 169 ± 30 

RV jitter HIRES (ms _1 ) d . . . 15.6 6.0 

Secondary eclipse parameters 

T s (BJD) 2455419.205 ± 0.058 2455200.132 ± 0.138 

T s 14 0.1289 ±0.0126 0.1424 ± 0.0107 

T s ,i2 0.0142 ± 0.0023 0.0183 ± 0.0074 

Planetary parameters 

M p (Mj) 0.626 ±0.037 0.778^;^ 

R p ( r j) l-212io!oM l-l 07 ioloia 

C(M P ,R P ) e 0.46 0.54 

p p (gem" 3 ) 0.44±0.09 0.71 ±0.18 

log 9p (cgs) 3.02 ±0.06 3.20 ± 0.07 

a (AU) 0.0434 ± 0.0007 0.0667 ± 0.0008 

T eq (K) 1384 ± 52 1260±45 

f 0.044 ± 0.003 0.077 ± 0.007 

F per (lO^rgs-^m^ 2 ) « .... 9 - 16 ±a92 6 - 90 ±i:oi 

F ap (10 8 ergs- 1 cm" 2 ) s .... 7.58 ± 1.34 4.72±°'g2 

(F) (10 8 ergs _1 cm -2 ) g S-^l^oe 5 - 69 ±o:75 

a T c : Reference epoch of mid transit that minimizes the correlation with the orbital period. 
BJD is calculated from UTC. T14: total transit duration, time between first to last contact; 
T12 — T34: ingress/egress time, time between first and second, or third and fourth contact, 
k Values for a quadratic law, adopted from the tabulations by Claret ( 2004) according to the 
spectroscopic (SME) parameters listed in Table [71 

G Lagrangian orbital parameters derived from the global modeling, and primarily determined 
by the RV data. 

d This jitter was added to the Keck/HIRES measurements only. The NOT/FIES observations 
are consistent with no jitter. 

c Correlation coefficient between the planetary mass M p and radius Ft. p . 

f The Safronov num ber is given by 6 = | (V cac /V orb ) 2 = (a / Rp) {M p / M 4 ) (see 
IHansen fc' Barman 2003). 

g Incoming flux per unit surface area, averaged over the orbit. 
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Fig. 9. — Mass-radius diagram of currently known TEPs. HATNet planets are show as green triangles and TEPs from other surveys 
are shown as blue circles. HAT-P-28b and HAT-P-29b are shown as large light blue triangles. The Solar system planets are shown as filled 
gray squares. The insert in the top left corner is a zoom of the region where HAT-P-28b and HAT-P-29b reside. Isodensity curves (in 
g cm -3 ) are plotted as dotted lines . Over laid are planetary 1.0 Gyr (brown, the upper set of lines) and 4.5 Gyr (green, the lower set of 
lines) isochrones from Fortncy ct al. (2007) for H/He-dominated planets with core mases of M c = Mgj (solid) and M c = 10 Mgj (dashed), 
interpolated to the solar equivalent semimajor axis of HAT-P-28b. Given the estimated age of 6-ltft Gyr, HAT-P-28b may be slightly 
larger than what is allowed by the theoretical models. HAT-P-29b is consistent with a H/He dominated planet with negligible core mass 
but also with a 1 Gyr model with a core mass of M c = 10 Mgj. 
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